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ABSTRACT 



We report on interferometric imaging of the CO J=l-0 and J=3-2 Une emis- 
sion from the controversial QSO/galaxy pair HE 0450-2958. The detected CO 
J =1-0 line emission is found associated with the disturbed companion galaxy not 
the luminous QSO, and implies Mga^Ha) ~ (1 - 2) x 10^° M©, which is > 30% of 
the dynamical mass in its CO-luminous region. Fueled by this large gas reservoir 
this galaxy is the site of an intense starburst with SFR ~ 37OM0yr~^, plac- 
ing it firmly on the upper gas-rich/star-forming end of Ultra Luminous Infrared 
Galaxies (ULIRGs, Lir > IO^^l©). This makes HE 0450-2958 the first case 
of extreme starburst and powerful QSO activity, intimately linked (triggered 
by a strong interaction) but not coincident. The lack of CO emission towards 
the QSO itself renews the controversy regarding its host galaxy by making a 
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gas-rich spiral (the typical host of Narrow Line Seyfert 1 AGNs) less likely. Fi- 
nally, given that HE 0450-2958 and similar IR-warm QSOs are considered typical 
ULIRG^ (optically bright QSO) transition candidates, our results raise the pos- 
sibility that some may simply be gas-rich/gas-poor (e.g. spiral/ elliptical) galaxy 
interactions which "activate" an optically bright unobscured QSO in the gas-poor 
galaxy, and a starburst in the gas-rich one. We argue that such interactions may 
have gone largely unnoticed even in the local Universe because the combination 
of tools necessary to disentagle the progenitors (high resolution and S/N optical 
and CO imaging) became available only recently. 

Subject headings: galaxies: active — galaxies: ISM — galaxies: starburst — 
ISM: molecules — quasars: individual (HE 0450-2958) — galaxies: mergers 

1. Introduction 

HE 0450-2958 is an optically bright (M^ = -25.8) IR-selected QSO at a redshift of 
z=0.286 (Low et al. 1988), and one of the few in the local Universe whose position in 
the (60;um, 25yum) versus (lOOyum, 60fim) far-IR color-color diagram lies in the area be- 
tween Ultra Luminous Infrared Galaxies (ULIRGs) and optically bright QSOs defined by 
-0.8 < a(100,60) < 0.4 and -2.0 < a(60,25) < -0.8 (Canalizo & Stockton 2001). Its 
complex environment, marked by tidal interactions, was identified early with ground-based 
(Hutchings, & Neff 1988) and /f^T imaging (Boyce et al. 1996). These efforts revealed 
a strongly interacting system consisting of the QSO and a disturbed galaxy ~ 1.5" away, 
setting this QSO/galaxy pair on par with typical ULIRGs in terms of morphology, and IR 
luminosity (Lir ~ 5 x 10-*^^ L0). 

The scarcity of such IR and optically luminous QSOs made HE 0450-2958 an early 
favorite candidate for undergoing a ULIRG—> (optically bright QSO) transition (Hutchings, 
& Neff 1988; Ganahzo Sz Stockton 2001) in a scenario first outlined by Sanders et al. (1988). 
The latter was proposed after FIR color-color diagrams have identified a population of dust- 
enshrouded QSOs via their AGN- heated dust component (de Grijp, Miley, & Lub 1987), 
and it was motivated by the similar bolometric luminosities and space densities of optically 
powerful QSOs and ULIRGs in the local Universe. It involves two gas-rich galaxies whose 
strong dynamical interaction induces a starburst and fuels an AGN within heavily dust- 
obscured environments. Then, as star formation uses and disperses their large molecular gas 
reservoirs, their initially cool Spectral Energy Distribution (SED) of dust emission changes 
towards one dominated by warm AGN-heated dust with much warmer IR colors. Eventually 
an unobscured optically bright QSO emerges out of the original ULIRG. The large dust and 
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molecular gas reservoirs found in optically selected QSOs where they fuel intense starbursts 
(e.g. AUoin et al. 1992; Haas et al. 2000; Evans et al. 2001) certainly support an AGN- 
starburst link, with dynamical interactions/mergers as the likely trigger for both activities, 
suggested also by optical studies (e.g. Canalizo & Stockton 2001). In this context the pivotal 
role of the ULIRG^QSO transition objects can then be understood as the "markers" of a 
special evolutionary stage during which the QSO emerges from a still ongoing dusty starburst, 
and lasting only a fraction of already short molecular gas consumption timescales (< 10^ yrs). 

HE 0450-2958 acquired further singular significance when Magain et al. (2005), after 
carefully subtracting the AGN emission from an ACS/HST image, failed to find the QSO 
host galaxy expected from well-established and tight correlations of (host galaxy)-quasar 
properties (e.g. McLure & Dunlop 2002; Floyd et al. 2004). This potentially pivotal result 
has been upheld by independent analysis (Kim et al. 2007), and follow-up VLT spectroscopic 
observations weakened the hypothesis for a significantly dust-enshrouded QSO host galaxy 
(Letawe et al. 2007). The possibility of an ejected (black-hole)-l-(fuel) system as responsible 
for a "naked" QSO, its serious consequences and observable signatures have been discussed in 
several papers. Summarizing their results here, such a configuration can be produced either 
by a Newtonian three-body "kick" exerted on a lighter black hole by a black hole binary 
residing in a gas-poor elliptical (Haehnelt, Davies, & Rees 2005; Hoffman & Loeb 2006), or by 
the recoil of a coalesced binary of spinning black holes due to an asymmetric gravitational 
wave emission (Haehnelt et al. 2005; Loeb 2007). The stage for both scenaria demands 
a strong galaxy interaction, and both have important consequences for gravitational wave 
detections by future detectors such as LISA, and the hierarchical build-up of supermassive 
black holes in galaxy centers (Haehnelt et al. 2005 and references therein). 

In the present work we report on the observations of CO J=l-0, the prime molecular 
line used to trace metal-rich H2 gas in galaxies, using the Austrahan Telescope Compact 
Array (ATCA), and the J=3-2 line with the Smithsonian Submillimeter Array (SMA). The 
paper is organized as follows: a) we present the observations and their analysis (section 2), 
b) estimate the molecular gas and dust mass of the system and discuss the properties of 
the companion galaxy (Section 3), c) present the implications regarding the AGN and its 
host (section 4), and d) discuss possible ramifications for the ULIRG— >QSO evolutionary 
scheme (section 5). Throughout this work we adopt a cosmology with Ho — 71 kms~^ Mpc~^, 
flu = 0.27 and Qa = 0.73, for which the luminosity distance of HE 0450-2958 at z=0.286 is 
Dl = 1460.4 Mpc and 1" 4.28 kpc. 
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2. Observations, data reduction, and imaging 
2.1. The ATCA observations 

We used the ATCA to image the CO J=l— hne emission (z^rest — 

115.2712 GHz) to- 
wards the QSO/galaxy system HE 0450-2958 during four periods in April, August 2006 and 
May, August 2007 utihzing the hybrid H214 and H168 configurations which combine an- 
tennas on both the east- west track and north track with maximum basehnes of 168 m and 
214 m respectively. These special configTirations provide good brightness sensitivity while 
achieving good u-v coverage within short (~6 hrs) observing sessions. This is important 
for mm observations at the ATCA site where atmospheric conditions usually limit useful 
observing (i.e. with low system temperatures and good phase stability) to time intervals 
significantly shorter than typical full synthesis sessions (~ 12 hrs). For the mean redshift 
z = 0.2864 of the pair (Canalizo & Stockton 2001), this CO transition is redshifted to 
i/(z) = 89.608 GHz, well within the tuning range of ATCA's 3 mm receivers (available in 5 of 
the 6 antennas). The correlator setup consisted of two IF modules of 64 x 2 MHz channels 
each, and frequency resolution of Az/jcs ~ 2.2Az/ch = 4.4 MHz. The effective velocity cover- 
age was AV(iFi+iF2) ~ (—300 +270) kms"^, after flagging bad edge channels and channel 
overlap (=7 channels) between the two IF modules. Pointing was checked hourly, and typical 
offsets were ~ 5" (rms), while the field of view of ATCA antennas at this frequency is 32" 
(HPBW). The array phase center was placed at the AGN's position, marked by its radio 
core at 8.6 GHz (the CI component in Feain et al. 2007). 

Dedicated wideband mm continuum observations were conducted during periods in May, 
August 2006, and August 2007 using the two IFs tuned in series around a center frequency 
~ 94.53 GHz (well away from the CO hne), and the same array configurations (two tracks 
with H 214, three tracks with H 168). For these observations the phase center was placed 3" 
away from the QSO as to avoid any (highly unlikely) small correlator DC offsets masking 
as a weak continuum source at the phase center. Typical system temperatures during line 
and continuum observations were Tsys ~ (200 — 350) K (including atmospheric absorption), 
estimated with hourly vane calibration performed per antenna. 

Amphtude/phase cahbration was obtained with 0454-234 observations interleaved with 

HE0450-2958 in lmin/(3-5)mins calibrator/source intervals, and passband calibration was 
achieved by observing 1921-293 or 2223-052. The absolute flux density scale was set by 
bootstrapping that of 0454-234 from Uranus and/or Mars observations in each track. These 
yielded S3mm(0454 — 234) = (1.9 ± 0.28) Jy, where the ~ 15% error represents the flux scale 
uncertainty of the data, obtained as the dispersion of S3min(0454 — 234) over the 2006, 2007 
observing periods (and thus is an upper limit since it contains also any real source variability). 
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2.2. The SMA observations 

A search for the redshifted CO J=3-2 line emission {urest = 345.796 GHz) was con- 
ducted using the 8-eIement SMA0 interferometer on 9th of December 2006 with a median 
atmospheric opacity at 225 GHz of T225 ~ 0.08. This transition is redshifted to 269.102 GHz, 
covered by the lower sideband (LSB) of the 345 GHz receivers placed at 268.892 GHz, while 
the upper sideband was offset by 10 GHz from this center frequency. The bandwidth covered 
by each is 2 GHz while the spectral resolution was 3.25 MHz for a total of 768 channels. 
This was smoothed so that a single CO J=3-2 map was created, averaged over a velocity 
range of 570 kms~^ centered at the redshift of the CO J=l-0 line emission. 

The SMA extended configuration was used, yielding a beam of ~ 1.20" x 0.94", PA=16.9° 
at the observed frequency and the u-v coverage attained. The mean system temperature was 
Tsys = 490 K (DSB), and the total on-source integration time ~ 5.07 hrs. Calibration of the 
antenna gain variations was done by observing the calibrators 0522-364 and 0455-462, at 
regular 10 minute intervals. Saturn and 3C454.3 were used to calibrate the bandpass, while 
Uranus was observed for flux calibration, with an expected uncertainty of ~20%. The MIR 
package was then used for the calibration of the visibility data. A search for 275 GHz 
continuum yielded an upper limit of S275GHZ < 8mJy/beam {3 a), over a ~ 3.936 GHz 
bandwidth (both sidebands), while the noise in the line map was arms ~ 6.6mJy/beam. 

2.3. Data reduction and imaging 

We used the in-situ ATCA phase monitor operating on a 230 m east-west baseline (Mid- 
delberg, Sault, & Kesteven 2006) to identify the periods when atmospheric phase noise at 
90 GHz over this baseline (comparable to the longest present in our datasets) exceeded ~ 30° 
(rms) and rejected the data taken during those periods. The remaining visibilities were then 
edited for any remaining large temporal amplitude/phase jumps. The final V(u, v) visibility 
dataset has a residual phase rms of cr^ ~ 20°, in accordance with a mean coherence factor 
Pcoh = (V(u, v))/S = e"'^'*'^'^'**^-'''^ ~ 0.95, estimated from imaging the visibilities of the cal- 
ibrator with an input flux density S (and antenna amplitude/phase solutions smoothed to 
Tavg ~ 2 X Tcycie(QSO ^ Calibrator) ~ 12 mins intervals). The u-v range of the final dataset 
is ~ (12 75) kA for both line and continuum. 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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Imaging was done using MIRIAD's task INVERT with natural weight used for max- 
imum point-source sensitivity in our line and continuum maps. CO J=l-0 emission was 
detected independently in both the II214 and H 168 array datasets, with peak flux densities 
differing by ~ 15%, consistent with the calibration uncertainties. An image of the CO J=l-0 
emission was produced by multi-frequency synthesi^ over the whole band, and combining 
both configurations. Deconvolution using the Clark algorithm was then performed. The 
resulting CO map and mm continuum image of the same region are shown in Figure 1, and 
the spectrum corresponding to the peak CO emission is shown in Figure 2. 



2.3.1. Source characteristics 

The velocity- averaged CO J=l-0 image yields the highest S/N possible from our data 
(Speak/o"rms ~15, Figure 1), and thus allows optimum determination of the position of the 
peak emission and the source size. A gaussian source model yields an excellent fit (i.e. after 
subtraction of the model, the noise in the residual image over the source area is arms = 
0.45 mJy/beam, as in the rest of the map area) where the peak CO brightness Speak = 
(6.7±0.45) mJy/beam is located at (Aa, A5) = (1.15"±0.12", -1.05"±0.10") relative to the 
AGN's radio core at the phase center. The small positional errors are those expected from the 
high S/N ~ 1/2 (Obcam)(S/N)-i ~ 0.09" relative to the phase center). The angular 

proximity of the amplitude/phase calibrator to HE0450-2958 (|Ak| ~ 6.5° away) allows 
also for excellent absolute astrometry by limiting 56'bas = (^B ■ Ak)/B ^ (50bas/27r)(0bcam) 
(the uncertainty due to the phase error 5(f)has ~ (27r/A)(5B • Ak)) to ~ 0.086" (for typical 
calibration errors of the baseline length of |5B| ~ 1mm). Thus we are confident that the 
peak CO emission is not located at the QSO position hut ~ 1.56" ± 0.17"(~ 6.7 kpc) away. 

The total flux density of Si*°*^(CO) = (10 ± 2) mJy emerges for a source for which 
the gaussian fit gives: 9maj = 3.95" ± 0.30" and = 2.80" ± 0.22" (PA=57° ± 9°). 

Deconvolving the restoring beam (3.21" x 2.14", PA=72°) yields an intrinsic size of 2.5" x 1.5" 
(PA=35°), larger than any expected "seeing" disk. The latter, deduced from images of the 
amplitude/phase cahbrator, is ~ 0.30" — 0.45" (consistent with the expected A6'seeing ~ 
[a0(rad)/27r] V8h^2 (Obeam) ~ 0.34", for a residual ~ 20° and (Obeam) = v^©!©^ ~ 2.62"). 

No continuum emission is detected at 94.5 GHz in either the starburst or the QSO 
position (S94GHZ < 1.4 mJy {3a), Figure 1). All relevant data, and the characteristics of the 
CO and dust emission in HE 0450-2958 are summarized in Table 1. 



^Fourier transform of the V(u, v) visibilities per channel, then averaging the resuhing channel maps 
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3. The molecular gas in HE 0450-2958 

The CO J=l-0 emission overlaid with the 8.4 GHz radio continuum is shown in Figure 2 
where a good correspondence with the region C2 associated with the companion galaxy 
(Feain et al. 2007) is evident. Regriding the observed CO emission and its deconvolved 
model, and then overlaying them onto the HST/ACS optical frame of HE 0450-2958 further 
underlines its association with the companion galaxy rather than the QSO (Figure 4). 

The corresponding H2 mass is given by M(H2) = XcoLco(i-o); where Lco(i-o) = 
/av /a TbdAdV is the velocity/area- integrated line brightness temperature at the source 
reference frame, and Xco (in M0(Kkms~^ pc^)~^ units) is the CO-H2 luminosity- mass con- 
version factor. Using standard derivations (e.g. Solomon et al. 1997) we have 



Leoa-o,^3.25xlO'(l..r.(^)-(^)^(te^X) Kk.s-pc^. (I, 

For the observed / Sco(i-o) dV (Table 1), Dl = 1460.4 Mpc and i/co,rest = 115.271 GHz (the 
hne rest-frame frequency), the latter gives Lco(i-o) — (2.3 ± 0.46) x 10^° Kkms~^ pc^. This 
is close to the maximum CO luminosity observed for starbursts, beyond which it remains 
constant across redshift (Frayer et al. 1999; Evans, Surace, & Mazzarella 2000), a possible 
indication of a self-regulating process occuring in extreme (maximal?) starbursts. Using 
Xco ~ 1 M0(Kkms~^ pc^)~^ (deduced for local ULIRGs, Downes & Solomon 1998), we 
obtain Mgai(H2) ~ 2.3 x 10^° M©, which places the companion galaxy in this galaxy/QSO 
pair firmly in the ULIRG category, at its most gas-rich upper end, and similar to starbursts 
found at high redshifts (e.g. Greve et al. 2005). 

A lower limit in the H2 gas mass can be obtained by assuming that the ^^CO J=l-0 
emission is optically thin. Using the derivation in Bryant & Scoville (1996) and assuming 
LTE, the new conversion factor becomes 



M(H2) 
Lco(i-o) 



0.0s 



J(TO - J(Tbg 
J(Tk) 



n -1 



M 



o 



Kkms~^ pc^ ' 



(2) 



where To = Ei/ke - 5.5 K, J(T) = To (e^^/^ - l) , Tbg = (1 + z)Tcmb ~ 3.5 K (CMB 
temperature at z=0.286), gi = 3 (degeneracy factor of n=l level), Z ~ 2(Tk/To) (partition 
function), and [CO/H2] ~ 10~^ (for solar metallicity) . For typical star forming gas where 
Tk ~ 40 - 60 K, the last equation yields (xS5o°^)Tk ~ 0.55 Mq (K kms"^ pc2)-\ and thus 
Mgai(H2)„,i„ ~ 1.25 X 10^0 Mq. 
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3.1. Constraints on the moleculeir gas excitation 

The CO J=3-2 observations yield a limit of Sco(3-2) < 3Nb^/'<5Sco(3-2)AV, where 
(5Sco(3-2) ~ 6.6mJy/beam is the noise of the CO J=3-2 map, Nb ~ 3.3 is the number 
of SMA beams corresponding to the CO J=l-0 emitting region and AV ~ 570kms~^. This 
gives Sco(3-2) ^ 20Jykms~^, and Lco(3-2) ^ 9 x 10^ Kkms~^ pc^, which corresponds to a 
(3-2)/(l-0) brightness temperature ratio of R32 < 0.39. This is rather low for molecular gas 
in starbursts where (R32) ~ 0.65 (e.g. Devereux et a. 1994), though within the range found 
in such galaxies (e.g. Yao et al. 2003). We examined an extensive grid of models with our 
Large Velocity Gradient (LVG) code (based on work by Richardson 1985), restricted by the 
upper limit on R32, and the condition Tkin > T^ust ^ 50 K (since the gas cools less efficiently 
than the dust, and turbulent /photoelectric heating heats the gas more). Most conditions 
reproducing this ratio have (n(H2)) ~ (10^ — 10^) cm~^ (volume-averaged H2 gas density), 
rather low for star forming H2 gas. This may in turn signify that not all of the CO (1-0)- 
emitting gas is involved in star formation (and thus some of it is not CO(3-2)-bright). Given 
the minimum R32(min) ~ 0.22 found in ULIRGs (Yao et al. 2003), deeper SMA observations 
of CO J=3-2 at a lower resolution (better matching the CO J=l-0 source size) should detect 
this hue in HE 0450-2958 and shed hght on its star-forming molecular gas phase. 



3.2. The companion galctxy: a typical ULIRG 

The mm/IR dust continuum SED of HE 0450-2958 from IRAS and our mm/sub-mm 
data (Figure 5) is typical for warm ULIRGs and IR-selected QSOs, where the presence of 
warm dust was the prime reason they were considered as ULIRG — > QSO transition objects 
in the evolutionary scenario linking these two classes. A cool/warm dust emission SED fit 
yields tJ""^^ = 45 - 55 K and T^™^ ~ 175 - 194 K (for emissivities of /5 = 1 - 2). Adopting 
p = 1.5 as our working value yields tJ;^°°'^ = 48 K, mJ;^°°'^ ~ 10^ Mq, and Lpm = 2.1 x lO^^ Lq 
(the luminosity of the cool dust component). A starburst in the gas-rich companion galaxy 
is a natural source of this large far-IR luminosity of HE 0450-2958, with a star formation 
rate of SFR ~ 1.76 x 1O-^°(Lir/L0) Moyr"^ ~ 37OM0yr-\ and a star formation efficiency 
of SEE ~ LFiR/Mgai(H2) ~ (90 - 165) (Lq/Mq), typical values for ULIRGs (e.g. Solomon 
et al. 1997). Moreover assuming most of the cool dust mass in this system residing in the 
companion galaxy, we obtain Mgai(H2)/M|j'|^g°''* ~ 125 — 230. These values arc well within 
the range found in IR-luminous galaxies (e.g. Sanders Scoville, & Soifer 1991, for the Xco 
adopted here), adding confidence to the estimated Mgai(H2) as a good measure of the bulk 
molecular gas mass present in HE 0450-2958. 

The results above are certainly in good accord with evidence of significant star formation 
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and reddening in the companion galaxy, and no such activity or dust obscuration evident 
in a putative QSO host galaxy (Magain et al. 2005; Letawe, Magain, & Courbin 2007). 
They contradict the recent view of Kim et al. 2007 regarding the nature of the companion 
galaxy, but also set HE 0450-2958 apart from typical ULIRGs with double nuclei. The latter 
have molecular gas mass ratios of ~ 1 : 1 — 2 : 1 (Sanders & Ishida 2004; Evans, Surace, 
& Mazzarela 2000), and their AGNs (when present) reside in the most gas- rich member of 
the interacting/merger pair (e.g. Evans et al. 1999, 2002). Subtraction of the CO emission 
model leaves no residual such emission in the location of the QSO with S^p^'^VSqso —host ~ 5 
(3(7), which corresponds to Mgai(H2)/MQso-host(H2) ^5:1 (for a common velocity range 
and Xco value) . A different (higher) Xco factor for the molecular gas in the QSO host could 
reduce the aforementioned asymmetry in H2 gas mass distribution, while some CO emission 
at the QSO's location at velocities outside our velocity coverage of ~ 570kms~^ could be 
missed. However the relatively small velocity range traced by the ULIRG, the QSO, and the 
highly excitated ionized gas around it (Mcrritt et al. 2006) make the latter possibility rather 
small. Sensitive wideband observations of e.g. CO J=3-2 with good brightness sensitivity 
would nevertheless be very valuable in the search for any QSO-related molecular gas. 



3.3. The dynamical mass of the companion galeixy 

The line profile towards the peak CO emission shows signs of a rotating disk (Figure 3), 
though the lower S/N per channel across the band makes this far from certain. The CO 
source model shows good overall agreement with the galaxy's optical image (Figure 4), and 
at z = 0.2864 where 2.5" x 1.5" — > 10.7 kpc x 6.4 kpc, its dimensions are comparable to the 
H2-rich parts of typical spiral disks (Regan et al. 2001; Heifer et al. 2003), including the 
Milky Way, but with ~ 10 — 20 times more H2 gas mass. 

For an underlying disk geometry the source angular dimensions correspond to an in- 
clination angle cos(i) = ^minor/^major = 0.6 =^ 1 ~ 53°. The enclosed dynamical mass then 
would be 



For R ~ Rmajor/2 ~ 5.35 kpc, VS'^ = (100 - 150) km s'^ (see Figure 2), it is Mdy„ ~ 
(1.9 — 4.4) X 10^° M0. Thus the galaxy's molecular gas mass amounts to > 30% of the 
dynamical mass enclosed in its CO-emitting region, making this ULIRG a very gas-rich 
system. 
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Inadequate resolution and/or low S/N makes it notoriously difficult to discern a gaseous 
disk from other possible configurations, and this has led to ambiguities regarding the presence 
of large gas-rich spirals, especially in the distant Universe (e.g. Genzel et al. 2003; Tacconi 
et al. 2006). In our case the resolution is inadequate to discern velocity/position patterns 
within the source and the underlying dynamical configuration could certainly be different 
(e.g. two orbiting compact gas reservoirs). Nevertheless for the intrinsic size and velocity 
width of our source most dynamical configurations would produce similar or even smaller 
Mdyn values (e.g. Bryant & Scoville 1996; Genzel et al. 2003), and thus the ULIRG in 
HE 0450-2958 would remain a very gas-rich object in terms of its Mgai(H2)/Mdyn ratio. 



4. Implications for the QSO's host galetxy and its AGN 

A suggestion for a Narrow Line Seyfert 1 (hereafter NLSyl) AGN (which have lower 
mass black holes), made by Merrit et al. (2006), could reconcile the lack of a massive 
spheroid host galaxy reported by Magain et al. (2005). However this rests on the weak 
premise of a good analogy between luminous QSOs and the much less luminous Seyfert 
galaxies, unsupported by the fact that NLSyl AGNs are typically hosted by gas-rich spirals 
with luminous starbursts, while optically luminous QSOs with My < —24 reside mostly in 
gas-poor massive elliptical hosts (Floyd et al. 2004). For example in I Zw 1 (the closest 
QSO and a prototypical NLSyl), a spiral disk hosts a massive cirumnuclear starburst at a 
radius of ~ 1.9 kpc (Schinnerer, Eckart, & Tacconi 1998; Staguhn et al. 2004) fueled by 
a large molecular gas reservoir. None of that is evident in the vicinity of the AGN in the 
HE 0450-2958 system whose bulk of the molecular gas and starburst activity are found in 
the companion galaxy instead. Indicatively, at z~ 0.286 our CO 1-0 observations could 
detect the molecular gas mass of the ~ 10 times less far-IR luminous I Zw 1 at a S/N~ 5. 
Finally, well-established empirical AGN- (host galaxy) correlations also predict a hitherto 
absent prominent elliptical (Magain et al. 2005), and thus this important issue remains open. 



4.1. The AGN-related molecular gas: could it be ejected? 

The two-component dust continuum SED fit (section 3.2) yields tS^^°'^(AGN) = 184 K, 
mJ^^^"^^ (AGN) ~ 5 X 10^ M0 and L^nid-m = 2.6 x 10^^ L© (the luminosity of the corresponding 
dust component), but unlike the typical scenario this AGN-heated dust reservoir now lies 
outside the ULIRG where the bulk of the dust and molecular gas mass resides. Using the 
same gas/dust ratio as in the companion galaxy gives Magn(H2) ~ (0.6 — 1.15) x 10'^ Mq. 
This amounts to the mass of a few GMCs, and such molecular gas quantities are found 
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within < 100 pc of AGNs in the local Universe (e.g. in the Sy2 galaxy NGC 1068, Planesas, 
Scoville, & Myers 1991; Schinnerer, Eckart, & Tacconi 1999). 

The disturbed nature of the HE 0450-2958 system and the extreme starburst in the gas- 
rich companion galaxy ~ 6.7 kpc away from the QSO argues in favor of a strong dynamical 
interaction. Could then a black hole, along with the aforementioned gas mass, have been 
ejected out of its host galaxy during such an event? In order to answer this and check the 
consistency of such proposals against the gas mass associated with this particular AGN, an 
estimate of its black hole mass is first needed. A lower limit can be derived from the fact 
that the dust within an AGN-bound gas reservoir will partly obscure it (giving rise to the 
distinct warm IR colors), and thus 

(AGN) _ 47rcGmH fcCEddMBH _ „ „ ^ „12 r / 



= =3.3.10" U [r^^ j L„ (4) 

where fc = l|^'~^'~~'''/Lagn < 1 is the fraction of the AGN's intrinsic luminosity intercepted 
by the dust and re-radiated at IR wavelengths, CEdd = LAON/LEdd (LEdd is the Eddington 
luminosity limit), (Tt is the Thomson cross section, and Mbh the black hole mass. Since 
LiR^^^ = 2.6 X 10^^ Lq, the last relation yields, 

MBH~0.8(fceEdd)"'lO'Mo. (5) 



For CEdd ~ 0.8 (typical for Mv(QSO) < -25, Floyd et al. 2004), and fc = (1 + Ropt.m)"^ 
0.71 (where Ropt,iR = (Lopt/LiR)AGN ~ [i'B^i^'^^\B)]/[uQQ^ji'^'-^^\miin\)] ~ 0.4, Ganalizo & 
Stockton 2001), it is Mbh ~ 1-4 x 10^ M©. This is consistent with the low values advocated 
by Merritt et al. (2006), but given that ours is only a lower limit, it does not settle the issue 
of the black hole mass (and thus of the expected host galaxy) raised by Magain et al. (2005). 

If such a black hole is the recoiling remnant of a coalesced black hole binary then, 
following Loeb 2007, it can carry a gas disk with a mass of 

M... ^ 3 . 10^„-«.»,-»« (^) ' ' M„ (6) 

where a ~ 0.1 is the disk viscosity parameter, rj = (eBH/0.1)(LAGN/LEdd) ~ 0.8, and Vcj 
is the ejection velocity of the recoiled merged BH product. For high ejection velocities of 
Vej > 1000 km s~^ (e.g Hoffman & Loeb 2006) Equation 6 yields the maximum M^isk that can 
be carried out by the recoiling black hole of Mdisk < 2 x 10^ [Mbh/(10^ Mq)]^'^ Mq. Given 
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that Mbh ^ 10^ M0, the latter comfortably encompasses the molecular gas associated with 
the AGN in HE 0450-2958. 



4.1.1. How long could the AGN remain luminous? 

Apart from partly obscuring the AGN in HE 0450-2958 via its dust content, the molec- 
ular gas mass in its vicinity can in principle "fuel" it for 



-fuel 



Magn(H2) 

Lagn/ (c^ cbh) 



~ 1.4x lO^eBHfc 





-1 


'Magn(H2)" 


1 10" Le ) 




IO^Mq 



yrs, 



(7) 



where cbh is the black hole fuelling efficiency. For a typical eBH ~ 0.1 and the estimated 
Magn(H2) the latter yields Tfuei ~ (2 — 4) x 10^ yrs. The shortness of this timescale, and the 
difficulty of fueling a black hole once it is outside its host galaxy, may point to fundamental 
limitations in the observability of such events, a plausible reason for their hitherto scarcity. 



4.2. The origin of the radio continuum 

The system's radio continuum contains low-luminosity jets emerging from a radio-quiet 
core located at the QSO position (Feain et al. 2007). An important question arising from 
that study is how could the location of HE 0450-2958 near the FIR/radio correlation be 
explained if the starburst occurs in the companion galaxy while most of the radio emission is 
of AGN origin. Clues for a possible answer may be offered by the case of 1821-1-643, another 
radio-quiet but IR and optically luminous QSO (My ~ —28, LpiR ~ 9 x 10^^ Lq) residing in 
a giant elliptical. This was the first radio-quiet QSO discovered to have classic FR I radio 
structure with low brightness diffuse jets (Papadopoulos et al. 1995), and a high brightness 
temperature compact core (Blundell et al. 1996), indicative of a black hole-based central 
engine. Despite the obvious AGN contribution in its radio continuum it did not deviate 
decisively from the far-lR/radio correlation until the discovery of a much larger, 300 kpc- 
sized, low-brightness FR 1 radio structure (Blundell & Rawlings 2001). A similar scenario 
may apply for HE 0450-2958, and the discovery of such a jet morphology in this system could 
also offer independent clues regarding the black hole. 

Interestingly, FR 1-typc radio structures may originate from a precessing jet axis due 
to binary black holes that have yet to coalesce (Blundell & Rawlings 2001). Such binaries, 
residing in gas-poor ellipticals (which allows for their long survival against gas-dynamical 
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friction), are a fundamental prerequisite in the scenario of ejected AGNs as the outcome of 
a three-body kick during gas- rich/gas-poor galaxy interactions (Hoffman & Loeb 2006). 

5. HE 0450-2958: not an ULIRG^QSO transition object, implications 

The hereby discovery of HE 0450-2958, considered an archetypal ULIRG— >QSO transi- 
tion object, as a strongly interacting pair of a gas/dust-poor galaxy (marked by a residing 
unobscured and optically powerful QSO) and a gas-rich extreme starburst raises the prospect 
of such interactions being misclassified in the context of the popular ULIRG^QSO evolution 
scenario. More such cases would then signify an important variation in the standard evo- 
lutionary picture in which starburst and QSO activity, while triggered by a common cause 
(a strong galaxy interaction/merger), do not necessarily involve two gas-rich progenitors 
with comparable amounts of molecular gas (Sanders & Ishida 2004). Such gas-poor/gas- rich 
galaxy interactions would not conform to the standard ULIRG— >QSO emergence scenario, 
though they can still ignite a starburst in the gas-rich progenitor (Di Mateo et al. 2007). 
They may also be the observational answer to the expectations of current galaxy-formation 
theories that predict the so called "wet-dry" mergers between gas-rich and gas-poor progen- 
itors to be the most common in the Universe (e.g. Springel et al. 2005). 

The demanding combination of observational tools it took to reveal these intriguing 
aspects of HE 0450-2958 may come to exemphfy why such interactions have been missed 
until now. Moreover, it could also resolve the conflicting views regarding the host galaxies 
of the more powerful optical QSOs (My < —24) in the local Universe. Unlike the low 
luminosity AGNs such as Seyferts where a consistent picture of AGNs in centers of gas-rich 
vigorously star-forming spirals has emerged (e.g. Schinnerer et al. 1998; Staghun et al. 
2004), detailed if^T studies by Dunlop et al. 2003 (see also Floyd et al. 2004) conclude that 
the more powerful AGNs reside in typical massive eUipticals (R^/^ light profiles, old stellar 
populations). On the other hand Scoville et al. 2003 (see also Bertram et al. 2007) found 
large amounts of molecular gas in such systems and concluded otherwise. 

It is now clear that the latter conclusions are premature. Indeed all current studies of the 
hosts of luminous AGN lack the combination of: a) high resolution optical imaging along with 
specialized deconvolution techniques (allowing sensitive probing of the host galaxy properties 
against the bright QSO), and b) interferometric CO maps with good resolution, high S/N, 
and the excellent astrometry that helped identify the type of merger in the HE 0450-2958 
system. Indicative of those limitations is PDS 456, the most luminous QSO in the local 
Universe (also a ULIRG/QSO composite in terms of its SED), where despite sensitive high 
resolution optical and CO imaging observations, the QSO host galaxy type and the location 
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of the molecular gas with respect to its AGN remain unclear (Yun et al. 2004). 

The short gas consumption timescales by a starburst, which are shorter still for optically 
powerful AGNs residing in gas-poor hosts, suggest that these twin "beacons" of gas-rich/gas- 
poor galaxy interactions may always be in close proximity. The advent of sensitive mm/sub- 
mm arrays, along with the already available high resolution optical imaging capabilities now 
allows their study, and the ULIRG— >QSO transition objects are excellent candidates. 

6. Conclusions 

We performed CO line emission imaging observations of the enigmatic QSO/galaxy pair 
HE 0450-2958 using the ATCA (CO J=l-0) and the SMA (CO J=3-2), detecting the J=l-0 
line and placing an upper limit on the J— 3-2 line, our conclusions are as follows 

1. The CO J=l-0 line emission is associated with the companion galaxy not the QSO, 
it corresponds to Mgai(H2) ~ (1 — 2) x 10^° M© and amounts to > 30% of the enclosed 
dynamical mass within ~ 8kpc. This large gas reservoir fuels extreme star formation at 
a rate of ~ 37OM0yr~^, efficiencies of LiR/Mgai(H2) ~ (90 — 165)Mo/L0, and places the 
companion galaxy on the upper gas-rich and star forming end of ULIRGs. This property, 
along with the disturbed nature of the system, solidifies the strongly interacting status of 
this controversial galaxy/QSO pair. 

2. No molecular gas is detected towards the QSO itself with MQso-host(H2)/Mgai(H2) < 1/5, 
consistent with the little dust extinction or star formation activity deduced for its putative 
host galaxy. This weakens the case for a host galaxy typical for Narrow Line Seyfert 1 AGNs 
since these are gas-rich spirals with luminous circumnuclear starbursts. HE 0450-2958 is thus 
the first known case of a gas- rich extreme starburst and a powerful QSO, intimately linked 
(via a strong interaction), but not coincident. 

3. Most of the radio continuum in this system is due to low-brightness AGN-driven jets 
rather than star formation activity. Deeper radio continuum observations are needed to 
reveal their morphology and extent, which in turn can yield more clues regarding its AGN. 

4. The early prominence of HE 0450-2958 as an IR-selected QSO that exemplifies the 
ULIRG^ (optically luminous QSO) transition in a popular evolutionary scenario, and its 
hereby revealed nature as a strong interaction between a gas-rich and a gas-poor, possibly 
elliptical, galaxy opens up the possibiUty that some of these so-called transition objects are 
such interacting pairs instead. Such gas-poor/gas- rich galaxy interactions can activate a 
starburst in the gas-rich progenitor, and an optically luminous AGN in the gas-poor one. 
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Their composite dust continuum would then be typical of IR-warm ULIRG^QSO transition 
objects. Sensitive optical and CO imaging at high resolution can now reveal and study such 
dynamical interactions, and these transition objects are excellent candidates. 
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Fig. 1. — Left: Map of CO J=l-0 emission (N A- weighted, CLEANed) averaged over 
570kms~"^ (see section 2.3), and contours: (—2,2,4,6,8,10,12,14) x iTrms- Right: The mm 
continuum emission at 94.53 GHz with contours of (—3, —2, 2, 3) x o"rms- The noise in both 
maps is Crms ~ 0.45 mJy/bcam, and the restoring beam is shown at the bottom left of the 
CO image (©beam = 3.21" x 2.14", PA = 72°). The phase center is marked by the cross, 
the diamond marks the AGN position, and the star marks the peak of the CO emission (see 
Table 1). For the mm continuum imaging the phase center was set 3" to the north of the 
AGN's position to avoid any DC correlator offsets "masking" as a spurious weak source. 
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Fig. 2. — The spectrum corresponding to the peak CO J=l-0 hne emission (Figure 1). The 
velocity range is centered on redshift of the QSO/ULIRG system shown on the upper left. 
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Fig. 3. — ATCA CO(l-O) greyscale image with ATCA 8.4GHz radio continuum contours 
overlaid. Contour levels start at 200/iJy beam~^ and increase in a geometric series with a 
common ratio of ^2. The CO(l-O) and the 8.4 GHz beams are shown in the bottom left 
corner of the image. The cross marks the peak emission of the C2 radio continuum emitting 
region associated with the companion galaxy (Feain et al. 2007). 
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Fig. 4. — Left: CO J=l-0 emission overlaid on the deconvolved HST/ACS image of the 
HE 0450-2958 system, with the QSO contribution subtracted (from Magain et al. 2005). 

Contours arc: (—2, 2, 5, 7, 9, 11, 13) x o"rms, with (Trms = 0.45 mJy/beam. The beam is shown 
at the bottom left (HPBW=3.21" x 2.14", PA = 72°). Right: The deconvolved CO source 
model (see section 2.3.1) overlaid on the HST image after a A6o = 0.6" shift aligning the 
optical and radio AGN positions. The good correspondence of the CO emission with the 
companion ULIRG is evident. In the optical image the QSO (marked by the single bright 
pixel near the phase center) is surrounded by an AGN-excited gas cloud. The cross size at the 
phase center (0", 0") denotes the (dominant) HST astrometric uncertainties (~ 0.5" — 1.0"), 
and the star symbol denotes the peak CO J=l-0 emission. 
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Fig. 5. — Two component fit of the mm/IR dust continuum of HE 0450-2958 (see section 
3.2). The upper Umits at mm wavelengths were obtained from this work (see Table 1), and 
the IRAS fluxes at 12 yum, 25/i.m, 60 /xm, and 100 yum are from Moshir et al. 1990. The 
emissivity laws adopted and the resulting temperatures from the 2-component dust emission 
fit are shown in the upper left (see also 3.2). 
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Table 1. Properties of the HE 0450-2958 system 



Property 


AGN 


Companion galaxy 


a(J2000.0) 


04^^ 52"^ 30M 


+1.15" ±0.12'"^ 


5(J2000.0) 


-29° 53' 35.0" 


-1.05" ±0.10"^ 




0.2863 


0.2865 


Speak(l-O) (mJy/beam)'= 


< 1.35 


6.7 ±0.45 


Speak(3-2) (mJy/beam)= 


< 20 


< 20 


/ Sco(i-o)dV (Jykms-^)'= 


< 0.77 


5.70 ± 1.15 


/Sco(3-2)dV (Jykms-i)^ 


< 11 


< 20 


Intrinsic CO (1-0) source size 




2.5" X 1.5" 


M(H2) (xlOi°M0)<^ 




1.25 - 2.3 


S94GHZ (mJy)'= 


< 1.4 


< 1.4 


S275GHZ (mJy)'= 


< 8 


< 8 


LiR (xIO^^Lq)- 


2.6 


2.1 



'^Position of the peak CO J=l-0 emission with respect to the 
AGN (see 2.3.1) 

•^Canahzo & Stockton 2001 

'^AU hmits are Scr, point source assumed for AGN-related hmits 

'^See section 3 

^from the dust continuum SED fit in Figure 6 



